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Abstract 

The minimum thermal conductance versus supercell size (d s ) is revealed in graphene and boron 
nitride superlattice with d s far below the phonon mean free path. The minimum value is reached at 
a constant ratio of d s /L ~ 5%, where L is the total length of the superlattice; thus the minimum 
point of d s depends on L. The phenomenon is attributed to the localization property and the 
number of confined modes in the superlattice. With the increase of d s , the localization of the 
confined mode is enhanced while the number of confined modes decreases, which directly results 
in the minimum thermal conductance. 
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Several recent experimental and theoretical groups have realized tunable electronic prop- 
erties from different graphene and hexagonal boron nitride (h-BN) compounds, flfl] With 
the further development of the synthesis technique, it will be quite possible for experimen- 
talists to produce the graphene and h-BN superlattice (CBNSL). The superlattice has been 
proposed as a promising thermoelectric material (for review see eg. Ref. 4). It can signif- 
icantly manipulate the thermal transport property while only introduces neglectable effect 
on electron transport. As one of the most important features of the thermal transport in 
superlattice, the minimum thermal conductivity has been discussed by lots of researchers in 
the superlattice of various components. 5Ml2l| In all existing studies, the minimum thermal 
conductivity is achieved in superlattice with supercell size d s at about the phonon mean 
free path. In 2000, Simkin and Mahan explained the minimum thermal conductivity as 
the result of the competition between the particle and wave nature of the phonon modes 
in the superlattice. 6J] The particle theory was applied for superlattice with d s larger than 
the phonon mean free path, while the wave theory was adopted in case of d s smaller than 
the phonon mean free path. These two different theories make opposite prediction on the 
behavior of thermal conductivity with the increase of d s , and a combination of them results 
in a minimum thermal conductivity for d s in the order of phonon mean free path. 

In this letter, we apply the nonequilibrium Green's function approach (NEGF) to inves- 
tigate the thermal conductance versus d s in the CBNSLs, where d s is far below the phonon 
mean free path. We report a new minimum thermal conductance in the ballistic superlattice 
at d s /L ~ 5%, where L is the total length of the superlattice. The value of the minimum 
point d s is proportional to the total length of the superlattice. It can not be explained by 
existing theory that predicts a constant minimum point d s around the phonon mean free 
path. We demonstrate that the minimum thermal conductance is related to various physical 
properties of the phonon confined mode (CM) in the superlattice. 

Figure. [T] shows the configuration of a CBNSL. The size of the supercell is d s = 2 in 
units of the smallest translational period of the graphene nanoribbon (do = 4.26 A). The 
total length of the superlattice is L = 10. Due to large phonon mean free path, [3, JjJ] 
the phonon-phonon scattering can be ignored for phonon transport in each graphene and 
h-BN sheet segments of the CBNSL; thus it is in the ballistic transport region. We refer 
to this type of CBNSL as ballistic superlattice. The thermal transport in the CBNSL can 
be investigated by the ballistic NEGF approach. 15] The phonon thermal conductance is 
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calculated by the Landauer formula: 

~dn(u,T) 



°ph = — / du JtjjjT[uj] 

Z7T 



dT 



(1) 



where T[u] is the transmission for phonon with frequency u, and n(u, T) is the Bose-Einstein 



distribution function. The force constant matrix is obtained from the GULP 



16] with the 



Tersoff+UOOP potential. 13, ll8] We have recently developed this set of potential for the 
graphene and h-BN sheet. It is efficient, stable and quite suitable for the study of heat 
transport. 

The thermal conductance in CBNSL versus d s is shown in Fig.[2]at 300 K and 1000 K. The 
value of the thermal conductance is reduced by the thermal conductance in pure graphene 
nanoribbon with same size (cr ). The conductivity k is related to the thermal conductance a 
through k = crL/S with S as the cross section area, so the ratio of the thermal conductivity 
k/kq will be exactly the same as u/uq shown here. Results of CBNSL with L from 60 to 
120 are shown in the figure. We have also observed similar phenomena in the CBNSL of 
different width and of zigzag edge configuration. A common feature in all curves is that 
the thermal conductance decreases with increasing d s in small d s region; then a minimum 
point is reached at d s /L w 5%; the thermal conductance increases with further increasing 
d s after the minimum point. The inset of panel (a) shows that molecular dynamics result 
also observes a minimum point around d s /L pa 5% for L = 120. The value of minimum 
point d s increases with increasing L. For L = 120, the minimum point is achieved at d s ~ 6, 
i.e., 2.6 nm. This length scale is about two or three orders smaller than the phonon mean 
free path in the graphene and h-BN, which is around hundreds nanometers. Actually, a 
similar minimum thermal conductance was also observed in the ballistic isotopic graphene 
nanoribbon superlattice with d s around 2 nm. [19|] It is not the one predicted previously by 
Simkin and Mahan, which should be in superlattice with d s around the phonon mean free 
path. [6j 

Figure. [3] shows the phonon spectrum of the graphene nanoribbon (left panel), the CBNSL 
(center panel), and the h-BN nanoribbon (right panel). Inset in each panel displays unit 
cell of the structure. The abundant difference in the phonon spectrum between graphene 
nanoribbon and h-BN nanoribbon leads to lots of CMs in the CBNSL. The CMs can be 
confined onto either C period or BN period, which is indicated by C or BN in the center 
panel of the figure. For instance, the C-type CMs of highest frequencies are confined in C 
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periods, corresponding to the in-plane optical phonon modes in graphene nanoribbon. In 
these C-type CMs, only C atoms vibrate while B and N atoms are silent. For the particular 
CBNSL in Fig. [TJ the C-type CMs can be confined in each of the five C periods; thus there 
are five CMs with similar frequency. If the CMs are confined perfectly, i.e., all B and N 
atoms do not vibrate at all, then there will be no information exchange among these five 
CMs. As a result, they can not transfer heat current and have no contribution to the thermal 
conductance of the CBNSL. However, the confinement in real materials is not so ideal and 
the B, N atoms can join in slight vibration, which is much weaker than the vibration of 
C atoms. In this situation, the CM can exchange information with its neighboring CMs, 
which enables it to transfer heat energy from one C period to its neighboring C periods. 
This kind of non-perfect CMs have contribution to thermal conductance. As pointed out 
by Ghanbari and Fasol in 1989, the interaction between two CMs through the superlattice 



layer decreases rapidly with increasing layer thickness. [20| Hence, the interaction between 
two CMs becomes weaker for larger d s , which enhances the confinement of these modes. We 
refer to this mechanism as mechanism A. As a result of the mechanism A, the CMs make 
few contribution to thermal transport with increasing d s . Another mechanism during the 
increase of d s with fixed L is that the number of the C/BN interface will decrease, leading 
to smaller number of CMs. This mechanism is referred to as mechanism B. As a result of 
mechanism B, the thermal conductance increases with increasing d s . The mechanism A is 
more important in the small d s region, which leads to the decrease of thermal conductance 
before the minimum point in Fig. |2j In the large d s region, the CMs are perfectly confined, 
so the mechanism B has dominant effect, leading to the increase of thermal conductance 
after the minimum point. The competition between these two mechanisms accounts for the 
minimum thermal conductance. 

To give a more quantitative characterization for the above argument, we investigate the 
localization properties of the CMs. The perfect CMs have better localization property than 



the non-perfect CMs 



re localization property can be well described by the inverse partic- 



22] The IPR for a phonon mode is defined through its normalized 



ipation ratio (IPR).[2ll. 
eigen vector u: 

N / 3 \ 2 

i=l \a=l / 

where N is the total number of atoms, and i = 1,2,3 are the three Cartesian coordinate 



components. The value of P^ 1 can be determined straightforwardly for two ideal special 
cases. P~ l = 1/N for the three acoustic translational phonon modes, where all atom 
has the same vibrational displacement u = l/yN. P^ 1 = 1 for phonon modes localized 
onto a single atom. A phonon mode with smaller IPR is better heat carrier. Fig. HJ (a) 
shows the localization property of the longitudinal acoustic modes in CBNSL of different 
d s . For convenience, we have plotted the iV\P _1 instead of P _1 in this figure. In whole 
Brillouin zone, the NP~ Y is about 1, indicating a good mobility of this acoustic mode. The 
acoustic modes are localized more strongly at the Brillouin zone edge. Considering the 
higher frequency of the acoustic modes at the Brillouin edge, this phenomenon is similar 
to acoustic properties of continuous media, where the evanescent decay rates are larger for 
modes of higher energy. (20] The localization property of the longitudinal acoustic modes is 
almost the same in superlattice with different d s . We have checked that this is also true 
for other two acoustic modes. Panels (b) and (c) show the localization property of C-type 
and BN-type CMs. Insets display the vibration morphology of the CMs. We calculate the 
NP~ l of these two CMs in case of various supercell length d s . Different from the acoustic 
modes, the NP~ l of CMs is sensitive to the value of d s , and increases with increasing d s . It 
manifests the enhancement of the localization property for the CMs in CBNSL with larger 
d s . As a result, the CMs make fewer contribution to thermal conductance for larger d s . This 
calculation demonstrates the origin of mechanism A. 

As the perfect CMs could not transfer thermal energy, an even more straightforward 
route to explain the minimum thermal conductance in Fig. [2] is to count the number of 
perfect CMs in the CBNSL of different supercell length d s . Fig. [5] shows the percentage of 
CMs in CBNSL of different length. A perfect CM is found if the IPR of this mode fulfills 
the criterion > Pq ■ Four different values are chosen for Pq . 0.05, 0.1, 0.2, and 0.3. 
Obviously, the percentage of CMs is smaller in case of larger -P _1 - However, a common 
feature in all four panels is that the percentage of CMs has a maximum value. The above 
mechanisms A and B are both incorporated in this figure. On the left of the maximum point, 
mechanism A is important. With the decrease of d s , some CMs turn to be able to pass heat 
energy to neighboring CMs, leading to smaller IPR value for these modes. As a result, the 
percentage of CMs decreases. On the right of the maximum point, the mechanism B is more 
important than mechanism A, since all CMs are already perfectly confined either onto C 
or BN period. With the increase of d s from the maximum point, the number of interface 
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decreases, which leads to the decrease of the CMs percentage. The competition between 
mechanisms A and B gives rise to the maximum percentage of CMs. This maximum point 
explains the minimum thermal conductance in Fig. |2J 

It should be noted, however, that the interface mode is another important type of localized 
mode, which exists in the SiGe superlattice. 20j The frequency of the interface mode locates 
in the middle of the energy gap between the phonon spectrum of the two components in 
the superlattice. A large value of energy gap is required for the existence of the interface 
mode. This condition is satisfied in the SiGe superlattice, as Si and Ge atoms have large 
mass difference leading to large energy gap. For the CBNSL studied in this work, the energy 
gap shown in Fig. [3] is not large enough to generate interface modes. Indeed, we do not see 
this type of modes in our calculation. 

In conclusion, we perform a NEGF study of the thermal conductance in the ballistic 
CBNSL. We find a minimum thermal conductance in CBNSL with supercell size d s /L m 5%. 
The value of the minimum point d s is far below the phonon mean free path in graphene or 
h-BN. This minimum thermal conductance is not the one predicted by previous works where 
d s is about the phonon mean free path. We demonstrate that the CMs in the CBNSL play 
an important role in this phenomenon. On the one hand, the CM becomes more perfectly 
confined for larger d s , leading to smaller thermal conductance. On the other hand, the 
number of the CM decreases with increasing d s , resulting in larger thermal conductance. 
The former mechanism dominates the small d s region, while the later mechanism controls 
the large d s region. The interplay between these two mechanisms explains the ballistic 
minimum thermal conductance. 

In this work, the d s of the minimum thermal conductance increases with increasing L, 
which reflects the ballistic nature of the phonon transport in the superlattice. The minimum 
thermal conductivity explained by Simkin and Mahan[6j occurs at d s around the phonon 
mean free path, which indicates the diffusive nature of phonon transport. An interesting 
open issue is to investigate the crossover from the ballistic minimum point to the diffusive 
minimum point. 
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FIG. 1: (Color online) The configuration of a CBNSL with total length L = 10 and supercell length 
d s = 2 unit cells. Heat transport across the superlattice from left to right. 
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FIG. 2: The thermal conductance versus supercell length d s at (a). 300 K and (b). 1000 K. <jq 
is the value of thermal conductance in pure graphene nanoribbon. Inset shows the results from 
molecular dynamics simulation at 300 K for L = 120. 
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FIG. 3: (Color online) Phonon spectrum in graphene nanoribbon (left), CBNSL (center), and h-BN 
nanoribbon (right). Insets are the corresponding unit cells. The symbols C/BN are with respect 
to the confined modes confined onto carbon/boron nitride periods. 




FIG. 4: (Color online) IPR value for different phonon modes in CBNSL: (a), longitudinal acoustic 
mode, (b). C-type confined mode, and (c). BN-type confined mode. 
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FIG. 5: (Color online) The percentage of confined modes versus supercell length d s in CBNSL with 
different length. The confined mode is found by its IPR P^ 1 > Pq 1 , with Pq 1 increases from 0.05 
to 0.3 in panels (a)-(d). 
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